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Abstract 

We describe an effective field theory for atomic lasers which reduces to the 
Jaynes-Cummings model in the non-relativistic, single mode limit. Our action 
describes a multi-mode system, with general polarizations and Lorentz invari- 
ance and can therefore be used in all contexts from the astrophysical to the 
laboratory. We show how to compute the effective action for this model and 
perform the calculation explicitly at the one loop level. Our model provides 
a way of analyzing a many-particle, two-state model with arbitrary boundary 
conditions. 
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I. INTRODUCTION 



The canonical model for laser physics is the Jaynes-Cummings model. It describes a 
single mode oscillator representing a coherent electromagnetic field, coupled to a two level 
reservoir of atoms Jl[. The Jaynes-Cummings model is defined in momentum space in terms 
of the photon creation and annihilation operators a and a) for a single momentum mode 
K = Q/c and a single, unspecified polarization of the electromagnetic field. The quantum 
mechanical Hamiltonian is given by 

H = a^ahQ + -huji2cr z + Tig(a + a + aV_), (1) 

where T1U12 is the energy difference between the atomic states. The creation and annihilation 
operators satisfy [a, a'] = 1 and the sigma matrices satisfy [<x + , er_] = a z . This model is the 
natural candidate for studying the fundamentals of the interaction between matter and 
radiation in a laser in a wide variety of situations, but it has several shortcomings and it is 
important to understand how these may be resolved in a reasonable fashion. 

In this paper we present a new model which retains the essential simplicity of the Jaynes- 
Cummings model, but which repairs some of its limitations. One of our principal aims is 
to write down a real-space Lagrangian formulation for a two-state system in which spectral 
content and polarizations are fully general: this should not only allow us to use the full 
machinery of modern field theory with all its attendant methodology, but also admit the 
solution of problems with general boundary conditions, in contact with many particle reser- 
voirs. The theory makes gauge symmetries and the space-time structure clearer and leaves 
us free to use well-established path integral or Green function methods for computing the 
effective action. Finally, but not least importantly, it also bridges a cultural gap between 
the worlds of field theory and laser physics. 

The Jaynes-Cummings model is an idealized description of laser phenomena. As a single- 
mode theory it can not address boundary conditions [@,|J or time-dependent interactions (4| 
since, by the uncertainty principle, a single mode must be completely delocalized in space 
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and time. These are features characteristic of the micro-maser and of non-linear media. 
Almost all of the work on lasers is phenomenological and couched in momentum space. Real 
space methods were pioneered by Graham and Haken , but the closest attempts at con- 
structing a microscopic description of the laser come from Korenman's |J use of Schwinger's 
action formulation H . A recent letter makes some progress with this approach for semicon- 



ductor lasers flip] . These papers also deal with effective theories however. Our paper is no 
different in this respect: we use an effective interaction and effective field variables. Indeed 
it would be inconceivable to attempt to write down a theory in which every optically-active 
electron and background charge were dealt with explicitly. Rather we pose the question: 
what are the relevant degrees of freedom for the laser at the energy scales of interest? These 
are clearly the averaged atomic properties and the magnitude of electromagnetic field. 

A disadvantage with Korenman's analysis is his use of non-relativistic field theory. Ko- 
renman begins with the Schrodinger equation coupled to reservoirs and seeks self-consistent 
solutions for decays rates and line widths. But radiative corrections to the non-relativistic 
theory are beset with problems: acausal loop diagrams, such as those used in constructing 
the effective action, vanish owing to the absence of anti-particles (negative energy states) 
in the non-relativistic theory. This makes the non-relativistic theory alien to field theo- 
rists who are used to the language of Feynman diagrams and Green functions and, in any 
case, one would expect a physical system described by the Schrodinger equation to arise 
naturally from a more general relativistic theory in the low energy limit. There is then 
the issue of non-renormalizability: Schrodinger scalar field theory is more divergent than 
relativistic scalar field theory, owing to the dimension of the field variables, and is specifi- 
cally non-renormalizable in 3 + 1 dimensions. It therefore makes more sense to begin with 
a relativistic theory, which is renormalizable, and consider the non-relativistic theory as an 
approximation to this full theory. In addition, we expect that a relativistic theory is nec- 
essary to study astrophysical situations, where the motion of atoms could be relativistic at 
high temperatures, even when the emitted radiation is of low energy. We begin therefore by 
introducing the action for a relativistic two-state model. 



II. THE ACTION AND ITS INTERPRETATION 



Consider a system of neutral atoms, containing optically-active electrons, which endow 
the atoms with a dipole moment. The electrons will not be explicit degrees of freedom in 
our model, rather their presence will be taken into account by the availability of transitions 
between the two atomic states. A neutral atoms is therefore represented as a two component 
real-scalar field; the two components represent the lower (unexcited) and upper (excited) 
levels of the atom. Each level has a different effective 'mass', in relativistic terminology 
m a = m + E a /c 2 , where m is the atomic mass and the potential-energy of the level is E a . 
In SI units, the action has the following form: 

S = J dV x \hi 2 c^d^ a ){dM + \mlc A 4> a <j) a + -^F^F^ + P^(0)F^| (2) 

where a = 1,2 and P^(0) is a polarization tensor which is to be specified below. Our 
conventions are such that the Minkowski metric tensor has the signature — h ++ and 
we use symbols da x to represent an n-dimensional infinitesimal spatial volume element on 
a spacelike hypersurface and dV x to represent the n + 1-dimensional infinitesimal spacetime 
volume which is canonically written da x dx°\/— detg. In other words, a x is a spatial volume 
and V x is a spacetime volume. 

The essential physics of this model is determined by the form of the dipole interaction 
tensor P M1/ (0). Given that the dynamical degrees of freedom are represented by real scalar 
fields, we have only two choices for this quantity, as we discuss below. The form for such a 
dipole term is unfamiliar in a relativistic theory, so we allow ourselves to be guided by the 
non-relativistic limit and require that this limit be consistent with known results, namely 
the non-relativistic analysis of Korenman || and in turn the Jaynes-Cummings model [|TJ. 
In particular, in the non-relativistic limit, one should obtain an expression for P^ u of the 
form used by Korenman: 

P,u - 7#V£Vfc, (3) 
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for some constant, off-diagonal matrix 7)™. It is evident that this is a dipole induced tran- 
sition from the form of the operators, ip* is a creation operator for the field and ip is a 
destruction operator, thus the off-diagonal operator creates an upper state and destroys a 
lower state, or vice-versa. Moreover, the components 7oi of this matrix will be proportional 
to the electric dipole moment of the atom. One relativistic generalization which reduces to 
eqn. (|), is 

Pp, = ihl^e ab <j)ado(J)b, (4) 
where 7^ is a constant, anti-symmetric tensor and e a b is the two-dimensional antisymmetric 



Levi-Civita symbol [11]. This form is intuitively appealing because it seems to be related to 



the relativistic inner product: 

(0a, <Pb) = ihc 2 J da x ^(<f)* a d <f) b - (<9 o O0b). (5) 

Unfortunately, this form for P^ u raises some questions concerning renormalizability (see 
section V). It is non-renormalizable in 3 + 1 and 2 + 1 dimensions: in particular we expect 
new, higher derivative interactions to be introduced at each order in perturbation theory. 
Although low energy predictions are still possible in such theories, we avoid this problem by 
introducing another interaction 

Pp, = \i^r b <p a <p b , (6) 

where 7^ is antisymmetric in \i and v and 

(oA 



—ab 



\ 10 / 



is now symmetric in a and b. In the non-relativistic limit, this interaction differs from that 
in eqn. (|]) only by a factor of ih/mc 2 (introduced by d t ). It has the advantage of being 
marginally renormalizable in 3 + 1 dimensions and super-renormalizable in 2 + 1 dimensions. 
This will be discussed in more detail in Section V, in which the one loop effective action is 
computed explicitly. 



We now sketch a simple derivation of the non-relativistic limit for our action, in order 
to give a physical interpretation to the constant matrix / y fiu . The next section contains a 
further justification of this method based on the field equations. The limiting procedure is 
unambiguous up to redefinitions of the origin for the arbitrary energy scale. The simplest 
procedure is to first observe that the real scalar field 4>(x) may be decomposed into 

=0(+)(x) (7) 

where (j)^(x) is the positive frequency part of the field and (j)(~\x) is the negative frequency 
part of the field and (j)^ + '{x) = (cj)(~'(x))*. We now rescale the fields by the atomic mass: 

^\x) = 4tL, 0(-)(*) = 4±L (8) 

V2mc 3 v2mc 3 
In addition, we note that the relativistic energy operator ihd t is related to the non-relativistic 
energy operator %hd t by a shift with respect to the rest energy of particles: 

ihd t = mc 2 + ihd t . (9) 

This is because the non-relativistic Hamiltonian does not include the rest energy of particles, 
its zero point begins just above the rest energy. 

Integrating the kinetic term by parts so that (<9 M </>) 2 — > </>(—□)</> and substituting eqn. 
(|) into eqns. (H) and (||) gives, 

+ J A7 x dt^jF^(i, + <J")„(<J + </>*)„. (10) 



Here we have dropped the Maxwell part of the action to avoid clutter, since it has no 
non-relativistic limit. If we use the fact that tp a (x) is composed of only positive plane- 
wave frequencies, it follows that terms involving ip 2 or (tp*) 2 vanish since they involve delta 
functions imposing a non-satisfiable condition on the energy 5(mc 2 + hu), where both m 
and uj are greater than zero. This assumption ceases to be true only if there is an explicit 
time-dependence in the action, indicating a non-equilibrium scenario, or if the mass of the 

6 



atoms goes to zero (in which case the NR limit is unphysical). In the next section we perform 
a transformation of the field equations which decouples the positive and negative frequency 
modes, justifying this procedure in a more conventional way. We are therefore left with 



Snr = hm J da x dt \ % — (^*(d t ip a ) ~ (d t ip* a )ip a ) ~ ^* a H a %l) a - 4 ^ ^(^6 + Mb) 



(11) 

where the differential operator H a is defined by 

and we have redefined the action by a sign in passing to a Euclideanized non-relativistic 
metric. It is now clear that, in the NR limit c — > 00, the final two terms in H a become 
negligable, leading to the field equation 

H a il> a (x) + J ^J F ^M*) = m?l>a(x), (13) 

which is the Schodinger equation of a particle of mass m moving in a constant potential of 
energy E a with a dipole interaction. The dipole interaction term is not negligeable since the 
constant 7^ is of order c 3 as we shall show below. 

The space-time components 7°* can now be related to physical electric dipole moments 
for linear media in the following manner. From classical electromagnetism we have that the 
dipole energy density is given by P • E, where P is the dielectric polarization and E is the 
electric field. The dielectric polarization is related to microscopic displacements of charge 
by 

P = eoXeE = — (er) x no. density of charges. (14) 

If we use the quantum number-density ip*ip here we see that the dipole energy density is 
given by 

p E = (er) • E. (15) 
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Our non-relativistic Lagrangian is also an energy density, thus comparing these in the rest 
frame of the charges, and using the fact that F 0% = —E l /c we have 

I—^f^ = -L^FkW = ~ -^WV = "(er • E)^, (16) 
allowing us to identify 

7 °* = mc 3 (er*). (17) 

Note that m is the mass of an atom and not the mass of the polarized charges. The spatial 
components are normally zero in the laboratory frame, but in relatively moving frames 
they may be determined by a suitable boost transformation. 

III. RELATIONSHIP TO THE JAYNES-CUMMINGS MODEL 

We now wish to show rigorously, making explicit the dimensionless parameters that 
have to be small in order for the approximation to work that, in the non-relativistic limit, 
our model describes a two level atomic system interacting via a dipole interaction with 
an electromagnetic field. For atoms interacting with a single radiation mode, the Jaynes- 
Cummings model emerges naturally. The Lagrangian for our model is 

- C = )fc\d^ a f + X -mYcg + ^-F^ + Y v e ah <p a <p b d^A u . (18) 

which gives for the equation of motion, 



□ + ~h- J & + TT^ab^d^ = 0. (19) 



To take the non-relativistic limit we define two fields IF2 



/ vn (p 1 % Ti 

= V - ~ (0a + ^a) 



m a c 



2 m a c z 



The rescaling is necessary in order that the non-relativistic wave-functions have the right 
dimensions, with standard inner product. Clearly, if a is real, then ip a = \* a - Moreover, 
these definitions imply: 
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\f2m~ a 



' (^a + Xa) (21) 



* - (¥) 7sb» w '» - * J (22) 

from which one can deduce: 

2 

W« + *.) = - X«) (23) 

This field redefinition reduces the action to one that is first order in time derivatives. Using 
standard Legendre transform theory, we obtain the part of the Hamiltonian density involving 
the scalar fields 



H(n a , <j) a ) = Ti a (j) a - C 



~ 2c 2 k 2 " a ' 2 

where 



tt^ + ^V(V0 a ) 2 + Iml^l + ^-FeM, (24) 



^ = C T7TT = c Ma = -icm/ — r— - Xa) (25) 

0(<Pa) V ^ 

and we have written 7 • F = 7^ F^ v for brevity. Replacing ix a and a by their definitions in 
terms of ip a and Xa as given above, we get the Hamiltonian density, 

H(l(>a, Xa) = m a C 2 i> a Xa + J^-(9i(^o + Xa)<9i(^a + Xa)) + e ab {^ a + Xa)(^Pb + Xb) (26) 



and 



£ = n a <p a - H(lT a , (f) a ) 



i^Xa^a - m a C tpaXa 

h 7 • F 

-(di(i> a + Xa)di(lp a + Xa)) ~ e ab (j> a + X a)(lpb + Xb) (27) 



4m a 4 A /m a m 6 

up to total derivatives. We obtain the equations of motion by varying with respect to ip a 
and Xa- 

iTl^a = V 2 (^a + Xa) + m a C 2 ^ a + ^ = (lp b + Xb) 

2m a 2c 2 ^/m a m b 

fl 2 €ab ry . p 

ih Xa = ~ V 2 (^a + Xa) - ^a^Xa ~ ~ . i^b + Xb) (28) 

2m a 2c 2 v /m a m b 
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and write the Hamiltonian as a four by four matrix in the 2-component space {a, b} crossed 
into the 2-component space {ip,x} : 

n 2 v 2 



H, 



ab 



$ab ( -TT— + m aC 2 ) + ~ ^ , = € ab 



2m n 



+ 



. / h 2 V 2 \ 7 -F 

Oab —— + „ 9 ; =£ab 



2m a J 2c 2 v /m a m b 



2c 2 ^m a m b 
O 



(29) 



where 





(l 


0^ 




( l) 


= 






1 = 








- 1 ) 




v -i oj 



(30) 
(31) 

So far these equations are exact. Note that the terms proportional to O are non- 
Hermition and couple positive and negative energy states. We want to perform a similarity 
transformation (non-unitary) that will remove the operator O that couples ip a and Xa- We 
use an operator of the form Up = e tA where A has no explicit time dependence. We will not 
be able to find the required A exactly, so we assume that a perturbative expansion exists in 
which A is small. In this case, <f>' = e tA <f> and <f>' = H'cp' which gives, 



H' = e iA He~ iA 

= H + i[A,H} + . . . 



(32) 



We will in fact need three independent expansion parameters. In addition to the usual 
non-relativistic expansion parameters A a for the two atomic states: 

h 2 v 2 



mlc 2 



(33) 



we will need the coupling expansion parameter: 

7 • F 



A, 



(34) 



mc 2 ^m 1 m2C 2 

where m = \{m\ + TO2). We will also assume that (mi — m 2 ) << m. Since, as argued above, 
7 • F oc mc 2 e(r • E), it follows that 
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A, 



e(r-E) 



mcr 



(35) 



Thus A3 << 1 requires the dipole energy in the electric field to be much smaller than the 
rest energy of the atom. 
We can now expand A: 



A = Ai + A 2 + 0(A 2 



where 0(A 2 ) refers to a product of any two of the small expansion parameters, and 



Al = -XaSabPO 



(36) 



A, 



-M ab f30 



(37) 
(38) 



In the above 



I 01 ) 



(39) 



It is easy to verify that i[A,H] to leading order in A exactly cancels the terms in H that are 
proportional to O. This decoupling of the ip and \ modes yields a Hamiltonian which can 
be written as a two by two matrix in the space {a,b} acting on the column vector 



^2, 



In particular (dropping the prime): 



H = Hq + H, t 



int 



(40) 



where the free part of the Hamiltonian is: 

/ 

H = 



^ + m lC 2 

2mi 1 



V 







2m2 



+ 77T.2C 



(41) 



This result is almost the same as Korenman's and differs only by the fact that the kinetic 
terms have different masses — a consequence of the fact that we have chosen to view the 
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shifted masses m a as fundamental. As seen below, this only gives a higher order correction 
which can be neglected in the non-relativistic limit. The interaction Hamiltonian is: 



/ 



int 





7-F 



7-F 

2c 2 v / mim,2 





\ 



/ 



(42) 



\ 2c 2 ^/mini2 

After some algebra the Hamiltonian can be written in the following form: 



H 



h 2 v 2 



2m 



+ mc 2 



+ 



7 • F 


^0 


n ( 


2c 2 v /mim 2 




2 

V 



1 
-1 

,^v 2 



2mmic 4 





\ 

/ 





-h 2 v 2 

2mm2C 4 J 



(43) 



where we have defined the energy difference u\2 = ( m i — rn^c 2 /%. 

The first three terms have a very natural physical interpretation: The first term is the 
free Hamiltonian for the "collective modes" of the atoms (the term proportional to the 
mass is just a shift in the energy and not relevant), while the next two describe the energy 
splitting and the corresponding dipole interaction with the electromagnetic field. Assuming 
that hu 12 << mc 2 , the last term is an order A correction to the second term. It is therefore 
higher order in the non-relativistic expansion and consistency demands that we neglect it. 

We drop the term in the Hamiltonian that corresponds to the collective modes and write 
the remaining piece as the sum of two terms: 



(44) 



where a z is the third Pauli matrix and a + and a_ are the usual raising and lowering operators 
for the atomic states: 



<7_ 



V 00 / 
v io, 



(45) 



(46) 



12 



We can now make contact with the Jaynes-Cummings model by assuming a single mode 
electric field, linearly polarized in the x-direction, as would be found in a high-Q cavity of 
volume V, for example. In terms of the standard harmonic oscillator creation and annihila- 



tion operators, the field can be written ||13||: 



E = x£ n (a + a ] ) sin Kz (47) 

where £q = [hn/e V]^ is the "electric field per photon" for an electric field of frequency 
Q. In the above, z points along the longitudinal axis of the cavity and K = Q/c is the 
magnitude of the corresponding wave number. 

The interaction Hamiltonian now takes the form 

H mt = - 3 T °!^L (<t + + a_)(a + at) (48) 
= — l xSn sin Kz(a, + cr )(a + a t ) (49) 

where ^ x is the component of 7o« in the direction of the electric field. 

We can drop terms proportional to o~_a and cr+a). These terms correspond to the si- 
multaneous lowering of an atom and absorption of a photon, and the simultaneous raising 
of an atom and production of a photon, and we expect them to be suppressed. We can see 
that this is the case by looking at the evolution of the operators in the Heisenberg picture. 
Writing 



<7±(t) 


= a ± (0)e ±iujt 


a{t) 


= a(0)e- int 


a\t) 


= a^(0)e int 



we find that a_a and a + a) are proportional to e ± *( w+fii ) t and the other two products are 
proportional to e ±l ( w_n )* . We are interested in a system that is tuned close to resonance 
Vt m lo and therefore, in the random phase approximation, terms proportional to e ±1 ^+ n ) t 
will average to zero because of the rapid oscillation of the phase. The final result has precisely 
the form of the interaction term for the Jaynes-Cummings model: 
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Hint = hg(aa + + a 1 a J) (50) 



and we identify the Rabi frequency in our model as: 



he 3 \Jm\mi 

This corresponds to the usual Rabi frequency P"3 



9 = - r^ /-= ( 51 ) 



^ = -^%^sin^ (52) 

on making the identification: 

lx = (ex) x ^m 1 m 2 c 3 (53) 

which is consistent with the identification for 7 made in the previous section (eqn. (|TTD), 
apart from terms of order (rn\ — m 2 )/m. 

IV. THE EFFECTIVE ACTION AND ITS INTERPRETATION 

Having established a connection to the Jaynes-Cummings model, we no longer need 
to refer to it and we can focus entirely on the relativistic case. Quantum corrections to 
the relativistic model may be computed using standard field theoretical prescriptions. The 
effective action is a particularly elegant way of generating such corrections. Although our 
theory is already an effective theory, this does not invalidate the procedure of looking for 
corrections due to correlations in our chosen field variables. If such corrections were already 
accounted for, they would simply renormalize away trivially in a renormalizable theory. In 
a non-renormalizable theory, this is not possible by definition and thus the corrections are 
always of interest. The calculation of loop corrections in non-renormalizable field theories 



has been shown to give accurate physics. See ref. [0] for a discussion of this. 

The effective action is a generating functional for one-particle irreducible quantum cor- 
rections. It describes completely the behaviour of the averaged field (or background field) 
after all quantum fluctuations have been taken into account. If one defines the functional 
integral for a field $ by 
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W[J] = -z In J dn{§\ exp + J dV J$| (54) 

where fr = c = /io = eo = l, and the c- number average field by 

$ = ($) = ^ (55) 

then the effective action is given by the Legendre transform of W[J] which displaces the 
explicit dependence on the source of fluctuations J, in favour of a dependence on the average 
field itself. 

r[$] = W[J] - JdV J$ (56) 

The resulting object is in all senses an action for the average field. In the remainder of the 
paper we make use of the background field method to compute the effective action. We 
begin by dividing the field into an average part and a fluctuating part for convenience: 

where <p a is the average field and <p a is the quantum field which replaces the total </> as the 



variable of integration in eqn (54). This division may now be used as a basis for generating 



a perturbation expansion for the effective action. Our model for the laser contains two fields 
4> a and A^. We shall assume that the average external field = 0, so that A^ = jj- may 
always be gauged to zero in all physical results, provided only that the systems lives in a 
box with a simple topology. It is nevertheless useful to keep this quantity non-zero when 
using r[0, A^] as a generating functional, since derivatives with respect to the vector field 
allow us to easily calculate physical quantities of interest. We expand the action around 
these background fields 



S 



a + (f a , A^ + A^ = Sdassf^, A^] + S 2 [(p a ,Aft] + S int \^ a , (p a , A^], (58) 



where S^ss is the term composed purely of background fields, S2 is quadratic in the quantum 
field variables and S^t is the remainder. The effective action is then given by the one-particle 
irreducible part of 
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= Sg a ,A,] + (S int }+t((S int ) 2 } + ... (59) 

We compute the effective action in two stages. First we consider ffuctuatations in the photon 
field leading to an intermediary effective action T^. These can be dealt with exactly and 
this leaves us with a result for the dynamics of the atomic system with all photon degrees 
of freedom eliminated. This is action could then be used to describe the situation in the 
micro-cavity maser where the measurable degree of freedom are the atomic states, and the 
effects of the photons are only felt indirectly. Secondly we consider fluctuations in the 
atomic degrees of freedom, such as one would expect in a gaseous or solid state laser. This 
gives us the full effective action T[0J. 

We begin by considering the radiation field with a Lorentz gauge fixing term added and 
associated Lagrange multiplier 1/a, 

S M = J dV x [\f^F^ + P^F, U + ^(V^) 2 } . (60) 

Ghost terms may be absorbed into the functional measure in view of the trivial nature of 
the gauge field contribution. 

The functional integral over A^ may be performed immediately since it is Gaussian. 
Integrating by parts and shifting the quantum gauge field (the field of integration) A^ — > 
A^ — 2d u P flu , one obtains without modification to the functional measure, 



'M 



J dV x |^ [-□ 5; + (l - ±) c^cr] A v -2j dV x ,(d»PnD„Ad P Pn} (61) 



The integral over the gauge field A^ is now a Gaussian and may be dealt with by standard 
results. This results only in a constant addition to the effective action which may be renor- 
malized away by a shift of the arbitrary zero point for the energy scale. The result is the 
one-loop correction 

r « [P^\ = const + 2 J dV x dV x , (d^)D va (d p P^) (62) 

where the free photon Green function is defined by the relation 
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V" = 8,?8(x,i!). (63) 
Using this general result we obtain the first stage effective action for atomic field a : 

= J dV x [-□ + m*] a + 2 J dVMx^W^ix, x')^ c {x')Ux')} (64) 

and 

V abc \x, x') = !ZlZCdS P D ua (x, x')). (65) 

where we have introduced the short hand notation 7^ = ^ u e ab . 

To generate the second stage effective action, we expand the atomic variables about 
a background or external field. The physical significance of this step is the presence of 
measurable averages for the atomic variables in our system. 

Expanding around free fields and dealing with the interaction term as an expansion of 
the exponentiated action, we get, 

= T A [0J - i In f duM (e lS ^) exp {tS^ } (66) 
= (S int ) + i ((S iQt ) 2 ) + ...0(f) (67) 

The Feynman (time-ordered) propagator is defined by 

(ip a (x)(p b (x')) = -iG ab S(x,x'). (68) 

For the renormalizable V vertex we now obtain the part of the effective action which 
is quadratic in the background fields. The outstanding terms do not contribute to the self 
energy and therefore to the decay rates of the atomic levels. 

{Sink) = ~i J dV x dV x ,V abcd (x, x')(^ a (x)0 d (x')G bc (x, x') + ^ a (x)^ c (x')G bd (x,x') 

+ (f) b (x)(f) c (x')G ad (x, x) + (f) b (x)(f) d (x')G ac (x, x')) + disconnected (69) 



The matrix 7^ is off-diagonal but symmetrical in a, b, so we may write 

V abcd = \v 1212 = F 2112 = F 2121 = F 1221 = -r-TidX D„ a (x,x>))\ . (70) 
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This last result gives us 

(S int > = -4^7^/ dV x j dV x ,(dX D ua (x,x')) 

x {0 1 (x)^ 1 (x')G 22 (x,x') + 2 (x)0 2 (x')Gn(x,x / )} . (71) 

We are interested in the effective coupling constants of the quantized theory, which may 
be defined through derivatives of the effective action. These provide us with information 
about the decay rates or lifetimes of the atomic levels and corrections to the Rabi-flopping 
frequency The momentum-space structure of these quantites also illustrate how photon 
energies are related to the interatomic spacings etc. Specifically we wish to compute the 
diagonal scalar self-energy £ aa , whose imaginary part gives an indication of the decay rates 
of the levels, 

S aa (x, x') = - ffffi, , (72) 
5(f) a {x)5(f) a (x') 

the interaction vertex (or generalized coupling constant) 

r (xx'x")—e h ggl^i (73) 
and the photon self-energy or polarization tensor 

IT < X J)= (74) 



V. ONE-LOOP RENORMALIZ ABILITY 

We now verify the one-loop renormalizability of our theory. For the remainder of the 
paper, we choose natural units in which H — c — eo = fJ>o — 1- General arguments indicate 
that renormalizablity is connected with power-counting, or the dimension of the coupling 
constant. In this scheme there is only one scale of dimensions. Length and time are com- 
pletely equivalent and mass is the inverse of length. A dimensional analysis of the action in 
these units leads to the conclusion that both the scalar and vector field a and in n + 1 
space-time dimensions has engineering dimension 
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1-n 



[0J = [AJ=L— (75) 

A renormalizable quantum field theory is one in which all the infinities accrued by the calcu- 
lational procedure can be defined away by reinterpreting the coupling constants appearing 
in the action. This is possible only if the infinite terms are of the same form as the original 
terms in the action which contain the coupling constants. In a non-renormalizable theory, 
it is not possible to absorb all infinities with a finite number of redefinitions. 

For the moment we will consider both the non-renormalizable and renormalizable inter- 
actions, in order to contrast them: 

n-l 
2 



P, v : [<f*\ = L- 

P,u: [-r\=L*? (70) 



We may consider the case of both two and three spatial dimensions, since a laser often has 
an axial symmetry which reduces its effective dimensionality. For the first of the interactions 
in eqn. ( |76|) 7^ has the dimensions of L in 3 + 1 dimensions and in 2 + 1 dimensions. In 
both cases the interaction is non-renormalizable. The second interaction is more successful. 
In 3 + 1 dimensions, 7^ is dimensionless which implies that the theory is strictly (also 
called marginally) renormalizable. In 2 + 1 dimensions, 7^" has the dimensions of L~^, 
which implies that the theory is super-renormalizable. 

Using the second interaction, the Jaynes-Cummings model can be represented as a renor- 
malizable field theory given by the action in eqns. (Q) and @. We note that, although 
renormalizablity is often regarded as a critereon for choosing between field theories, it is not 
an infallible guide to their physicality. Quantum corrections to non-renormalizable theories 



are known to give accurate results in a number of cases [nfl. Moreover, we have a natural 
energy cut-off for the kinetic motion of atoms, namely kT. Our primary reason for choosing 
the renormalizable interaction is that it is easier to calculate quantum corrections in this 
case; the lack of an explicit time-derivative preserves Lorentz covariance. 
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A. Scalar self-energy 



We seek to calculate the one-loop self energy E(p) and vertex function T^(p,p',q) and 
show that these terms have infinite pieces that have the same form as the original interaction, 
and thus can be reabsorbed into the coupling constants, the masses, and rescaling factors. 
We use cutoff regularization since we ultimately want to use our model to study laser physics, 
which will involve the imposition of boundary conditions. We define the subtraction scheme 
by expanding around the mass shell. We expand in where 5 = ml — m\ and M 2 = 
\(m\ + m 2 2 ). 

The bare theory gives a propagator of the form, 

i 

iGaaip) 



P 2 + mlo ~ it 

which has a pole at p 2 = —m 2 a0 . We calculate the polarization tensor and use the Dyson 
equation to obtain a propagator of the form, 

iGaa(p) 



P 2 + m lo + Soa(p) -ie 
We define Am 2 = m 2 a — m 2 aQ and write, 

iG aa (p) 



p 2 + m 2 a + [E M (p) - Am 2 a } - it 
We choose 

Am 2 a = Ka(-m 2 a ) (77) 

so that the pole occurs at p 2 = —m 2 a which we call the physical mass. Thus, the propagator 
can be written, 

iG aa (p) = 



p 2 + m 2 a + [£ aa (j0 2 ) - Eaa(-n»a)] _ ie 

^aa(p) is divergent and a requirement for renormalizabilty is that we can write (after regu- 
larization) 
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Z aa (p 2 ) - ^ aa (-ml) = {p 2 + ml)f(A) (78) 
so that the propagator becomes 

iGaaip) = 2 J Z t . ; V = 1 + /( A ) ( 79 ) 
p 2 + m z a — %e Va 

These redefinitions are equivalent to the statement that we can add counterterms to the 
Lagrangian of the form, 

and absorb the infinities from the self energy S aa (p) in the mass shift Am 2 and the wave- 
function scaling factor Z<f, a . 

In a similar way, the interaction part of the Lagrangian C int gives rise to a bare vertex 
of the form, 



r<?> = 2e abV7o ^ 



where q is the incoming photon momentum and 7^ is the bare coupling constant. The one 
loop contribution to this vertex is divergent. We isolate the divergent part by performing a 
subtraction at the mass shell, 



^ ) (p,p',q) = T^ 



ras 



(1) 



where the subscript ms means that the external momenta are evaluated on the mass shell, 
and T is finite. A requirement of renormalizability is that we can write, after regularization, 

r ? ) L = r ? ) ( z r 1 - 1 ) ( 8 °) 

which means that we can absorb the infinite part of the one loop vertex graph into a redefi- 
nition of the coupling constant. This redefinition is equivalent to adding to the Lagrangian 
a counterterm of the form 
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In this section, we will calculate the one loop self energy and the one loop vertex function 
and use J77|), (|T5D, ( |79| ) and (§(]), to determine Am^, Z^ a and Zj. We start from the following 
expression for the self-energy of the scalar field 4>\. 



(2vr) 4 (P - ze) ((Jfe + p) 2 + m| - ie) 
where we have used the Feynman gauge a = 1 for the internal photon propagator. The 
self-energy for the field 02 will depend on mi in the same way. We rewrite the denominator 
using the usual Feynman parameter formula, 

dx- 



k 2 — ie (k + p) 2 + m 2 — ie Jo (k 2 (l — x) + [(k + p) 2 + m 2 }x — ie) 2 
We complete the square in the denominator and shift the integration variable k = I — px to 
obtain, 

dH r 1 l T l x + x 2 p T p x 



4*7 T/ V. 77^7 dx 



(2tt) 4 Jo (I 2 + a 2 - ie) 2 
where a 2 = m\x + p 2 x(l — x) and we have dropped the terms linear in I which give zero 
by symmetric integration. We do a Wick rotation so that the integration contour lies along 
the imaginary axis and make the change of variable, Iq = iU to obtain the Eucledian space 
integral, 

v a x f dH f 1 l T l\ + x 2 p T p x . . 

Sl1 = 47 7 "7 J^Y Jo dX (l 2 + a 2 ) 2 (81) 

The integral is infinite and we use cut-off regularization to render it finite. After regulariza- 

tion we can switch the order of integration and perform the / integration first. We consider 

the two pieces separately. First we evaluate the term proportional to l T l\ and call it Ef x . 

Under the integral sign we can replace l T l\ by \g T \l 2 (by symmetric integration) which gives, 

dH I 2 



(2tt) 4 (I 2 + a 2 ) 2 
where r y 2 = 7^7^ '. Doing the I integration gives, 

^ n = -f- [ 1 dx(A 2 + a 2 -2a 2 \n^] 
{Ztt) 6 Jo a 
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We expand in b = 5/M 2 and take only the leading order term. We calculate m 2 ) and 

G° 2 ) ~~ ^11 ( — m i)- The result is, 

<-v 2 2 A 2 M 2 

ZU-rnl) = t^(A 2 - -M 2 lnA, - — ) 
nK x> (2tt) 3V 3 M 2 9 

SL(P 2 ) - E^-m?) = -^[|(P 2 + m?) + \(P 2 + m?)ln^] (82) 

Next we have to calculate the term proportional to From ( [81]) we have, 

ci 4 / 1 



(2tt) 4 Z 2 + a 2 

where 7 2 A = 7^7^ ■ Doing the / integration we obtain, 

.,2 ,! A 2 



which gives, 



«)=|^v[i.4-i] 

S ;; ( p 2 ) - Eft-m?) = _^ P V[^i] (83) 

Thus, from fl77|), (fF|), (^|), fl82|) and fl83|) we obtain, 

. 2 7 2 „ 2r A 2 2, A 2 1 ^7 2 pVri, A 2 1.. 
Am = — '— M 2 — - -In— - - + 4 . . -In—? - - ( 84 ) 
a (2vr) 3 L M 2 3 M 2 9 7 2 M 2 L 3 M 2 9 JJ V ; 

Z * (27F [ T8 + 3M^ + 3^ 2 M 2_J (85) 

and Zfo = Zfa. Keeping only the divergent terms, Z^ 1 ^ is a wavefunction renormalization 
factor of the usual form, and the first two terms in Am 2 give an infinite shift in the mass 
term in the standard way. The fourth term in the expression for Am 2 corresponds to a new 
interaction in the Lagrangian at the one loop level of the form, 
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B. Interaction vertex 



Next we obtain the vertex renormalization constant from the one loop vertex correction 
shown in Fig XX. We obtain, 

r (i) = o,- y* bX //« f 

" 1x1 1 ^ J {2iry{k* -it){{p> -k) 2 + ml-it){{p-ky + ml-it) 

We rewrite the integral in terms of two Feynman parameters by using the expression, 

1 



with 



ABC Jo Jo a (xyA + x(l -y)B + (1 - x)Cf 

A = (p - k) 2 + mj - it 
B = (p - kf + m\-ie 
C = k 2 - ie. 
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Putting the external scalars on the mass shell we obtain, 
where 

b 2 = x 2 m\ + 5x 2 (l -y) + x 2 y(l - y)q 2 
We shift the integration variable 

I = k — (pxy + p'x(l — y)) 

and drop the terms linear in I which give zero by symmetric integration. We perform a Wick 
rotation so that the integration contour lies along the imaginary axis, and make the change 
of variable ko = ik\. The result is, 



a/3 



(2tt) 4 7o Jo * (P + b 2 Y 
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where, 

M a/3 = (p a xy + p' a x(l - y))(ppxy + p' p x(l - y)) 

We consider separately the terms proportional to l a lp and M a $. We will first do the 
integral for the term containing l a lp and call it T^ 1 . This term is divergent, and we use 
cutoff regularization. By symmetric integration we can write l a lp = \g a ^ 2 under the integral 
sign. Switching the order of integration and performing the / integration gives, 

r (i)/ = _4 7 W l x dx £ dy(ln± - |). 

We set q 2 = and do the integrals over x and y to obtain, 

We calculate the term proportional to M a p in the same way. Including this result, we obtain 
from O and (g§), 

x 7 2 / A 2 1 2 7 2 /3 W/ 3 (pV + pV) , q a g^\ 
Zl =1 -j2^[} ln Ip-2 ] + — [ ^- 2M 2 + 3i^ ] J (8?) 

The terms in square brackets represent contributions from new interactions of the general 
form, 

This higher-derivative term could become important in the strong field limit and in non- 
perturbative regimes. 

C. Photon polarization 

Finally, we consider the photon polarization tensor. We have, 
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We separate the denominators using the Feynman parameter technique and perform a Wick 
rotation as before. The result is 

where 

M 2 = m\{l -x) + m\x + g 2 x(l - x) 

Expanding around the mass shell, we isolate the divergent piece by setting q 2 = 0, which is 
equivalent to taking the first term in the expansion. The result is, 

n„(g)=4g a g^V^[ln^ + l] 
which leads to an induced interaction of the form, 



D. Renormalizability Revisited 

In the introduction it was claimed that a dimensional coupling constant matrix 7^ 
ensured the renormalizability of the model. However, the above calculations show that the 
one loop divergences require counter terms of the form 7 2 i ,(<9 M </>)(<9 !/ 0) and |7 AM 7 T ^i 7 A jU F Tiy . 
These terms may be thought of as multiplicative modifications to the scalar field kinetic term 
in the Lagrangian and to the value of /j,q and eo in the Maxwell part. They arise because 
the orientation of the dipole 7^ breaks the rotational invariance of the theory, which is 
then reflected in the quantum corrections. We are not obliged to add counterterms of the 
form 7 aA 7 / ^7 Miy 0a(<9a0fc)(<9/3 , 9/i^) since these new interactions yield finite results (at least to 
one-loop), but the appearence of such terms nevertheless indicates that they are an integral 
part of the structure of the relativistic theory and should therefore be considered too. 

The question then arises: is the theory, as given, renormalizable or not? We point out 
that, in a renormalized field theory, it is the renormalized values of the parameters which 
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are to be identified with the physical constants in an experiment. In our case 7^ is to be 
indentified with the dipole moment of an atomic system. In fact, the one loop divergences 
simply tell us that there are additional, relativistically covariant terms that are second order 
in derivatives of the fields that we could have added to the classical action. These terms 
correspond to relative permittivities and permeabilities. These are the only such terms 
which need to be added, and with the addition of these terms, the model would indeed be 
fully renormalizable. However, note that these two terms involve 7 2 , which we assume is 
small, so according to the assumptions on which we base our perturbative expansion they 
are probably negligable. Their physical significance is not redundant however: in the limit 
of large electromagnetic fields, very high kinetic energies and strong dipole couplings, these 
extra terms become significant and predict new physics to be identified with experiments. 
For the present paper, we take the pragmatic approach however and assume that such terms 
will not contribute signicantly. In effect we are renormalizing the new 7 2 couplings to zero. 
This is a significant improvement over the non-renormalizable choice of coupling, in which 
new, higher derivative interactions would appear at all orders in perturbation theory. In 
renormalization group philosophy, one would say that we are expanding our theory in a 
region of Lagrangian-space which is closer to a renormalization group fixed point. 



VI. CONCLUSIONS 

We have presented a relativistic model for the interaction of a two state atom with an 
electromagnetic field and verified that it reduces to the Jaynes-Cummings model in the 
appropriate limits. We have also shown how to compute higher order, quantum corrections 
and verified that the model is one-loop renormalizable. By identifying the renormalized value 
of 7 M i, with observed dipole moments, or the Rabi flopping frequency of known systems, we 
have a prescription for gauging the magnitude of corrections which lead to the onset of new 
physics. The decay rates of the atomic levels may be identified with imaginary contributions 
to the self-energy S aa , for which we are able to calculate an explicit expression, rather than 
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merely a formal expression as in Korenman's work. In the present paper, we have been 
mostly concerned with the self-consistency of our proposed model and have presented only 
a zero-temperature expression for the self-energy. In future work we shall compute the finite 
temperature self-energy, where the natural cut off Ac 2 is of the order of kT and obtain a 
more accurate gauge of the decay rate by looking at retarded (causal) boundary conditions, 
rather than the Feynman boundary conditions used here. It will also be natural to look at 
non-equilibrium systems, and extend our ananlysis to non-linear phenomena where some of 
the assumptions made in this paper begin to falter. 

Most laser systems are well described by non-relativistic physics. We consider the most 
important result of our paper to be the identification of a model which can be straightfor- 
wardly solved in real-space, with arbitrary boundary conditions, as well as in many-particle 
theories at finite temperature and non-equilibrium. The use of relativistic field theory sim- 
plifies calculations greatly compared to direct non-relativistic formulations. It also addresses 
quantum corrections at the level of the Lamb shift, where corrections are measurably sig- 
nificant in atomic systems |T5|Ji|] , and removes some of the arbitrariness of previous work on 



lasers by tying laser physics to a model which can easily be be analysed within the framework 
of a renormalization group philosophy. This is significant because it indicates which results 
are independent of the specific details of microscopic theory one chooses to work with. 

Our paper opens a doorway to the study of the statistical mechanics of photons and atoms 
in cavities and free space, a topic which we intend to pursue in later work. Interesting studies 
include the use of our model to study the micromaser with proper finite boundary conditions 
and partially reflecting surfaces, and in an expanding or contracting spherical cavity, as a 
toy model for light generation by bubbles in sonoluminescence, and porous silcon. 
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